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1 Introduction 

A wide class of the mathematical physics problems is connected with forced 
nonlinear Schrodinger equation (NLSE). Usually authors consider either the 
small amplitude solutions [3 E] or the asymptotics for perturbed soliton so- 
lutions ISlinilZ!- The problems on soliton generation usually are not discussed 
or solved numerically |H1 Ej • 

In this paper we show a mechanism of the soliton generation in nonlinear 
Schrodinger equation due to a small fast oscillating driving force: 

ia0^ + %^ + |^|2^ = £Ve^^^'', 0<e<l. (1) 

We investigate formal asymptotic solutions with a small amplitude of the 
order of O(e^). The leading-order term is the solution of linear nonhomo- 
geneous equation. We consider the one phase asymptotic solution with the 
phase function S/e^. Other terms of the order of 0{e'^) do not play any role 
in investigated phenomenon. 

The specific of the problem consists in existence of a local resonance. The 
resonance curve is defined by equation: 

deS+{d^Sf = Q, (2) 

where S is the phase function of the driving force. The amplitude of forced 
oscillations increases from 0{e'^) up to 0{e) in the neighborhood of this 
curve. 
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The increase of the amphtude is similar to a phenomenon of a local res- 
onance in the linear nonhomogeneous ordinary differential equation of the 
second order with the constant coefficients pH]- The amplitude increases 
linearly. The domain of this increase usually is called a resonance layer. The 
additional term of the order of the layer width will appear in the asymptotic 
solution after passage of the resonance layer. The order of residual 

part of the solution does not change. 

After resonance the leading-order term of the asymptotic solution has the 
order of 0{e). The spatial variable is defined by amplitude of the driving 
force. The relation between the order of the solution amplitude and typical 
scale of the independent variables is such that the leading-order term of the 
asymptotics is defined from NLSE. So specific relations usually appear due 
to specific asymptotic representation of the solution as was shown in well 
known works on asymptotic transfer to NLSE |3], ^Hj-PS]- In this work we 
show that such kind relation can be obtained due to resonance with driving 
force. 

Thus after resonance layer the leading-order term of the asymptotics has 
the order 0{e) and is defined from the Cauchy problem. The initial data 
are equal to the value of the amplitude / of driving force on the resonance 
curve. The further evolution of the leading term of the asymptotics does not 
depend on driving force. 

The effect of the solution increase after local resonance is well known 
in the ordinary differential equations ^U] and partial differential equations 
[TTl[T2j. Also it is well known the effect of appearance of NLSE in asymptotic 
constructions for the problems with strong dispersion, cubic nonlinearity and 
the appropriate relation between scales. This problem contains both of these 
effects. The passage through the local resonance leads to the increase of the 
amplitude. The modulation of the amplitude of the driving force and scales 
of the variables lead to NLSE behind the resonance curve. The similar results 
can be obtained for other problems on the resonance passage by driving force 
in systems with the cubic nonlinearity and strong dispersion. 

The considered problem is the model problem for the wide class problems 
on the local resonance passage by driving force. In order to show essence of 
this phenomenon we assume that the amplitude / of perturbation is a smooth 
function with respect to variable and phase function is = {e^OY /2. 

The result of this work is the following proposition 
Proposition. One phase asymptotic solution of the order of 0{e'^) in the 
domain 9 ^ —e^^ will have the order of 0{e) in the domain e^^ <^ 9 < 
Ke~^, K = const > 0. The leading-order term of the asymptotics 
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is determined from the Cauchy problem for NLSE 

.0 I ,2 

I Ut + Urcx +\ U \ U= U, 

u \t=o = (1 - i)y/nf{x). 



In particular, if the amplitude of the perturbation force has the form 

1 2?7exp(-z/i(a; + ?/o)) 



fix) 



'l — i)y/7r cosh(2?7(a; — xo)) 



then the leading-order term u of postresonance expansion is a soliton of 
NLSE, where xo,yo,fi,ri are soliton parameters. 

2 Forced oscillations 

In this section we construct the formal asymptotic solution for equation 
in preresonance domain. 

Lemma 1. In the domain 9 -C —e~^ the formal asymptotic solution for 
equation 0) with respect to base 0{e^) has a form 



e''A{ei,e''e)+e^A{ei,e^e)+ (3) 



exp(zS'/£:^), e ^ 0, 



the coefficients A of the asymptotics are determined from algebraic equations 

Equation is written in terms of fast variables ^ 9 but further for 
convenience we also use slow variables x = and t = e^O. 

The substitution Q into equation ((T)) gives a recurrent sequence of equa- 
tions for coefficients of asymptotics: 

- S' A= /, 

-S' A= -i At -d^x A, (4) 

,3 2 2 1 n 1 

-S' A= -i At -dxx A-\A\^ A. 

The solutions of this systems have singularities at t = 0. The constructed 
asymptotic solution can be represented as t ^ in the form: 



^{ei,t,e) 



p2 { l\ \ p^f^I_ _ III 
^t)^ \t^ t^ 
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f^'^fxx fxxxx l/P/ 3/ 



exp(iS'/£:^ 



The structure of singularities as t ^ —0 of coefficients in asymptotic rep- 
resentation (jni) allows us to determine the domain of validity for external 
asymptotics: 

Lemma 1 is proved. 



3 Internal solution 

In this section we construct the formal asymptotic solution for equation 
into resonance layer. The asymptotics of these type are usually called by 
internal asymptotics TB^. 

Lemma 2. In the domain — e"^ <^ 9 <^ e"^ the formal asymptotic solu- 
tion for equation (QJj with respect to base 0{e'^) has a form 

<^{i,T,e) = ew{ei,T)+e^w{ei,T)+e^w{ei,T) e ^ 0, (5) 

where t = eO . The coefficients w of the asymptotics are determined from 
ordinary differential equations (EP; ^Bil- 

The structure of singularities for coefficients in (jS)) as t ^ —0 allows us 
to determine internal scaled variable t = t/e = eO. 

The internal asymptotic solution is constructed in form The be- 
haviour of coefficients can be obtained by matching: 

w (e^, r) ~ 



exp{ir^/2}, r — oo. (7) 

exp{irV2}, r -oo. (8) 
The leading-order term of (0) is determined from equation: 

i /exp(irV2). (9) 
The solution of this equation can be written out in terms of Fresnel integral: 

w= -if {si) r exp{i0y2)de. (10) 

J —oo 
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exp{ir^/2}, r — > — oo. 



(6) 



w {ei, t) 



fxx ^ifx 



\f?f 



w {ei, t) ~ 



fx 



The first order correction w is determined by equation: 

•1 |0|90„ 

I + Wxx +1 W I W= 0, (11) 

the second order correction is determined by equation: 

I + +2\ w y w + w ^ w = {]. (12) 

In order to determine the behaviour of the asymptotic solution after res- 
onance passage we need to know the asymptotics of the coefficients of (jSJ as 
r — s> +00. Direct calculations lead us to formulas 



exp(2r^/2) 

IT 



where C\ = (1 — i)\/7i; 



^ (ee, r) = TW M)+ ^ oH) + qM)^^^^^^^ + (^(r- 



where 

W 1 = Ci/a;^: + |Ci/|^Ci/, 

il;oK)= lim f r [^..(£e,^)+|^(£^,^^)r^(£e,^)]o!^-(ci/xx+|ci/pci/)V 

T— +00 Y J —OO J 

gi{s^) = kifxx + ^2!/^/, fci and ^2 are constants; 
^ (ee, r) = ^ 2K) + r w M)+ ^ oH) + 92^ ^''^^'^'^^^ + 0{t~'), 

IT 

where 

2 / . // 1 \ r^i ,9 1 9 1 ^\ 
w 2{ei) = I \ {wi)xx + 2| Wo I wi + w qW ^ \ , 

w i{ei) = I \ {w o)xx + 2| Wo I Wo +2 w w _i wi - 
—2 w_iw Q Wi +2 w w _i w ^+ w Q w Q I , 

Wo= lim (i r [w xx{ei,e) + 2\w {ei,e)\^ w {ei,e)+ w ]de 

\ J —00 



— r W2 — r Wi 

and the function (72 (^^0 is smooth and expressed by f{e^). 

The asymptotic behaviour of the coefficients in (jSj) gives us the domain 
of validity for this solution: 

|r|<£:"^ or |6'| < 

Lemma 2 is proved. 
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4 Postresonance solution 

In this section we construct the formal asymptotic solution for equation 
in postresonance domain. 

Lemma 3. In the domain e^^ 6 < Ke^^ the formal asymptotic solu- 
tion for equation (OJ) with respect to base 0{e'^) has a form 

^(^, e,e)=eu {ei, e^O) + e\ u {e^, e^e)+ B e^O) exp{iS/e^)) + 
e% u (ee, s^e)+ Bi ^^0) exp{iS/e^)+ B-i e^O) exp{-iS/e^)) + (13) 
+e\ Bi {ei,e^d)exY>{iS/e^)+ B2 {ei^e^O) exY>{2iS / e^) + 

+ B-1 {ei,e^e)exp{-iS/e^)), 

m 

the coefficients Bk of the asymptotics are determined from algebraic equations 
\i4\ l and the coefficients u are solutions of the Cauchy problems for either 
NLSE or linearized NLSE U^-^. 

The substitution of (fT^ in equation ((T)) gives: 

e(^ - S' B -f^ exp{iS/e'^) +e'^(^i Ut + u^^ +\ u \^ 

Uut + u^^ +2\uf u + u'^ u + 

+i Bt exp{iS/e^)+ Bxx exp{iS/e^) + 
+2| M 1^ S exp{iS/e'^)+ u'^ B* exp{-iS/e'^) - 
-S' Bi exp{iS/e^) + S' |„i exp{-iS/e^)j 

I Z Mi + M^a; +2| M p U + M ^ U + 2| U p M + U ^ M * + 2 M I B p + 

+ M * i? ^ exp(22S'/e^) + 2 m * u5 axpiiSje^) + 
+2uu* B exp{iS/£'^) + 2 uuB * exp{-iS/£'^) + 
+2| u |2 Bi exp{iS/e^) + 2| w ^ B-i exp{-tS/e^) + 

+ u'^ B*i exp{-iS/e'^)+ u'^ Bli expiiS/e"^) + 

2 2 
+i Bi texp{iS/e^)+ Bi xx Q'x.p{iS / e^) - 

2 2 ^ 

-i B^i f exp(-iS'/£: )+ B^i xx^y-vi-^'S/ e ) + 

-S' Bi exp{iS/e^) + S' B-i exp{-iS/e^) + 2S' B2 exp{2iS/e^)) = 0{e'') 
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1 2 

Functions B, B± are determined from algebraic equation: 

- S' B= /, 



S' Bi= 'I Bt + Bxx +1 u \' B, 
S' B-i= -u^ B* 

r^/ 3 .2 2 , ,2 1 ^ ^ 1 1 Oil 2 2 ^ ,^ 

S' Bi= t Bi t+ Bi XX + \ u \^ B +2 u * uB + uu * B + u ^ B (14) 

,3 2 2 oil o 2 

B-i= i B-i t- B-i XX- uuB *~u^ Bl, 

,3 1 o 

2S' B2=u * B ' 



The matching with internal solutions give us the recurrent system of the 

12 

Cauchy problems for coefficients u, u, u: 

.0 I ,2 

Z Mf + Uxx +1 I U, 

u \t=Q = ci/(x); (15) 



2^4 + U^r^ +2 M -U + M M =0 



M \t=o =Wo (x); (16) 



■2 2 0,2 2 022^ i^i^O ,1,2 I2O 

t Ut + u^x +2\ u \ u + u u = — \ B \ u — \ u \ u — u u 
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u |i=o =w^o {x). (17) 

i 

Coefficients ij^ have singularities at t = 0. The domain of validity for 
solution (|T3|) is determined by inequalities: 

t > £ or 61 > 

Lemma 3 is proved. 

In order to determine the behaviour of solution ()13j) for large t > it is 

necessary to know the structure of the leading-order term u. It is easy to see 
that the behaviour of u is determined by amplitude of the perturbation force 
in original equation (pj). In particular, if the function f{eS,) has the form: 

,1 .X 1 2r7exp(-i/i(£^ + |/o)) 

J^^^' = 71 tT^^TTT — FT^' ^^l^ = const, 

(l-2)v^ cosh(2r7(e{-{o)) 

then to determine the leading-order term of asymptotics ()13|) it is necessary 
to solve the Cauchy problem for NLSE with initial data of the soliton type. 

We are grateful to L.A. Kalyakin and B.I. Suleimanov for stimulating 
discussions and for help in improving of the mathematical presentation of 
the results. 
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